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Abstract. We apply the statefinder diagnostic to the torsion cosmology, in which 
an accounting for the accelerated universe is considered in term of a Riemann-Cartan 
geometry: dynamic scalar torsion. We find that there are some typical characteristic 
of the evolution of statefinder parameters for the torsion cosmology that can be 
distinguished from the other cosmological models. Furthermore, we also show that 
statefinder diagnostic has a direct bearing on the critical points. The statefinder 
diagnostic divides the torsion parameter ai into differential ranges, which is in keeping 
with the requirement of dynamical analysis. In addition, we fit the scalar torsion model 
to ESSENCE supernovae data and give the best fit values of the model parameters. 
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1. Introduction 

The current observations, sucli as SNela (Supernovae type la), CMB (Cosmic Microwave 
Background) and large scale structure, converge on the fact that a spatially homogeneous 
and gravitationally repulsive energy component, referred as dark energy, accounts for 
about 70 % of the energy density of universe. Some heuristic models that roughly 
describe the observable consequences of dark energy were proposed in recent years, 
a number of them stemming from a certain physics [1] and the others being purely 
phenomenological [2] . Dark energy can even behave as a phantom and effectively violate 
the weak energy condition [3]. 

In various cosmological models, fundamental quantities are either geometrical 
(if they are constructed from a spacetime geometry directly) or physical (if they 
depend upon physical fields). Physical quantities are certainly mo del- dependent, 
while geometrical quantites are more universal. About thirty years ago, the bouncing 
cosmological model with torsion was suggested in Ref . [1] , but the torsion was imagined 
as playing role only at high densities in the early universe. Goenner et al. made a general 
survey of the torsion cosmology [5], in which the equations for all the PGT (Poincare 
Gauge Theory of gravity) cases were discussed although they only solved in detail a few 
particular cases. Recently some authors have begun to study torsion as a possible reason 
of the accelerating universe [Ej. Nester and collaborators [7j consider an accounting for 
the accelerated universe in term of a Riemann-Cartan geometry: dynamic scalar torsion. 
They explore the possibility that the dynamic PGT connection, reflecting the nature 
of dynamic PGT torsion, provides the accelerating force. With the usual assumptions 
of homogeneity and isotropy in cosmology and specific cases of the suitable parameters 
and initial conditions, they find that torsion field could play a role of dark energy. 

One of the motivation was to avoid singularities in the initial investigations of 
torsion cosmology [8]. However, it soon was found that non-linear torsion effects were 
more likely to produce stronger singularities [9] . The non- linear effects turn out to play 
a key role for the outstanding present day mystery: the accelerated universe. In the 
various PGT, the connection dynamics decomposed into six modes with certain spin 
and parity: 2^, 1^, 0^. Some investigations showed that 0^ may well be the only 
acceptable dynamic PGT torsion modes [10]. The pseudoscalar mode 0~ is naturally 
driven by the intrinsic spin of elementary fermions, therefore it naturally interacts with 
such sources. Consequently, it is generally thought that axial torsion must be small and 
have small effects at the late time of cosmological evolution. This is a major reason 
why one does not focus on this mode at the late time. On the other hand, the scalar 
mode O"*" does not interact in any direct obvious fashion with any known type of matter 
[TT] . therefore one can imagine it as having significant magnitude and yet not being 
conspicuously noticed. Furthermore, there is a critical non-zero value for the affine 
scalar curvature since O"*" mode can interact indirectly through the non-linear equations. 
The homogeneity and isotropy of cosmology have received strong confirmation from 
modern observations, which greatly restrict the possible types of non-vanishing fields. 
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Under the assumption of homogeneity and isotropy, 0"*" mode has only a time component 
and it can be specified as the gradient of a time-dependent function. Therefore, the 
cosmological models with the scalar mode offer a situation where dynamic torsion may 
lead to observable effect at late time. We emphasize again that one does not focus on 
the early universe, where one could indeed expect large effects (though their signature 
would have to be separated from other large effects), and substitutionally asks about 
traces of torsion effects at the late time of cosmological evolution [7]. 

Obviously, the fine-tuning problem is one of the most important issues for the 
torsion cosmology And a good model should limit the fine-tuning as much as 
possible. The dynamical attractor of the cosmological system has been employed to 
make the later time behaviors of the model insensitive to the initial condition of the 
field and thus alleviates the fine-tuning problem [12]. Furthermore, Nester et al [7] have 
shown that the Hubble parameter and a have an oscillatory form for the scalar torsion 
cosmology. 

The traditional geometrical parameters, i.e., the Hubble parameter H = a /a 
and the deceleration parameter q = —aa/d?, are two elegant choices to describe the 
expansion state of universe but they can not distinguish various accelerating mechanism 
uniquely, because a quite number of models may just correspond to the same current 
values of H and q. However, Sahni, Saini, Starobinsky and Alam [13] have introduced the 
statefinder pair {r, s}: r=d /aH^, s = (r — l)/3(g — 1/2). It is obviously a natural next 
step beyond H and q. Fortunately, as is shown in the literatures [TU [15] , the statefinder 
parameters which are also geometrical diagnostics, are able to differentiate a series of 
cosmological models successfully. Using the discussion of statefinder parameters in the 
scalar torsion cosmology, we explain easily why the present field equations modify the 
expansion of the universe only at late time. If the evolving trajectory of statefinder have 
a decelerating phase (g > 0) at early time, then we can understand why the expansion 
of the universe until z ~ few remains unchanged in the scalar torsion models. 

In this paper, we apply the statefinder diagnostics to the torsion cosmology. 
We find that there are some characteristics of statefinder parameters for the torsion 
cosmology that can be distinguished from the other cosmological models. The 
statefinder diagnostics show that the universe naturally has an accelerating expansion 
at low redshifts (late time) and a decelerating expansion at high redshifts (early time). 
Therefore, scalar torsion cosmology can avoid some of the problems which occur in 
other models. Especially, the effect of torsion can make the expansion rate oscillate 
when torsion parameter Oi > or ai < —1. Whether the universe has properties 
which are easier to explain within the scalar torsion context is a remarkable possibility 
demanding further exploration. The oscillatory feature of Hubble parameter had earlier 
been reported for the braneworld cosmology [16] and the quasi-steady state cosmology 
[T7] . We show that statefinder diagnostic has a direct bearing on the critical points 
of the dynamical system. One of the most interesting characteristic of the trajectories 
is that there are loop and curves with the shape of tadpole in the case of the torsion 
parameter — | < ai < 0. In this case, we fit the scalar torsion model to current type la 
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supernova data and find it is consistent with the observations. Furthermore, we analyze 
prehminarily the relevance for realistic observation of the found statefinder parameters. 

2. The equations of motion 



PGT [18] has been regarded as an interesting alternative to general relativity because 
of its gauge structure and geometric properties. PGT based on a Riemann-Cartan 
geometry, allows for dynamic torsion in addition to curvature. The affine connection of 
the Riemann-Cartan geometry is 

r^/ = r^/ + ^(T^/ + + T%) (2.1) 

where F^^, is the Levi-Civita connection and T^^ is the torsion tensor. Meantime, the 
Ricci curvature and scalar curvature can be written as 

+ ^{T^„^T^% + 2T,,<,T'^%) , (2.2) 

R =R + 2V^T^ + ^(T^.^T^'^'^ + 2T^,^T'^'"' - AT^T^) , (2.3) 

where Rf^i, and R are the Riemannian Ricci curvature and scalar curvature, respectively, 
and V is the covariant derivative with the Levi-Civita connection (For a detailed 
discussion see Ref. [7].). Theoretical analysis of PGT led us to consider tendentiously 
dynamic "scalar mode". In this case, the restricted expression of the torsion can be 
written as [7] 

2 

Tfj,,yn = -T[^gv]K, (2.4) 

where the vector T,, = T^,,, is the trace of the torsion. Then, the Ricci curvature and 
scalar curvature can be expressed as 

R^.u = R^u + ^(2V.T^ + g^.JWkT'') + ^(T^T, - g^.TJ"^) (2.5) 

R =R + 2V,,r^ - ^T^T^ . (2.6) 

The gravitational Lagrangian density for the scalar mode is 

L= -^R+^R^ + ^(T.^^T'^'^^ + 2T,.^T'«^^ - AT.T'^) , (2.7) 

where R = R'^'^ and ai is a torsion parameter. Consider that the parameter b is 
associated with quadratic scalar curvature term R^, so that b should be positive [7]. 
Therefore, the field equations of the scalar mode are 

V,R+^iR+^)T, = 0, (2.8) 
ao{R^, - ^g^^jR) = -{T^u + %u), (2.9) 
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where 7^,^ is the source energy-momentum tensor and T^^, is the contribution of the 
scalar torsion mode to the effective total energy-momentum tensor: 

6 y 

- - ^g^,R) . (2.10) 

Since current observations favor a flat universe, we will work in the spatially flat 
Robertson- Walker metric ds'^ = —dt'^ + a'^{t)[dr'^ + r'^{d6'^ + sin^ 6d(j)'^)], where a{t) is the 
scalar factor. This engenders the Riemannian Ricci curvature and scalar curvature: 

Rt' = 3- = 3{H + H^) , (2.11) 
a 

R/ = r/ = r 4> = ^ + 2^ = H + 3H^ , (2.12) 

• 2 

:r = 6(- + ^) = 6(if + 21/2), (2.13) 

where a{t) is the scalar factor, and H is the Hubble parameter. The torsion should 
also be only time dependent, therefore one can let Tt{t) = $(t) ($(i) is the torsion field) 
and the spatial parts vanish. The corresponding equations of motion in the matter- 
dominated era are as follows 

H = J^R-J^-2H\ (2.14) 
6ai 6ai 

$ = _^i?__^_3/7$ + l$2 (2.15) 

2ai 2ai 3 ' ^ ^ 

where n = ai — Oq and the energy density of matter component 

P = ^(^ + f){3H - $)^ - Ai?^ - 3a,H^ . (2.17) 

One can scale the variables and the parameters as 

t Ip'^H^^t, H IpH^H, $ ^IHq^, R I^HqR, 

«o ^pOo, ai -> /pfli, /i /J/i, b Ip'^Hq%, (2.18) 

where Hq is the present value of Hubble parameter and Ip = VSttG is the Planck length. 
Under the transform (12.181) . Eqs. fl2.14p - fl2.16l) remain unchanged. After transform, 
new variables t, H, $ and R, and new parameters oq, ai, fx and b are all dimensionless. 
Obviously, the Newtonian limit requires ao = —1. 

For the case of scalar torsion mode, the effective energy-momentum tensor can be 
represented as 

=-3;.i7^+A(i?+^)(3i7-$)^- (2.19) 

f/ = T/ = f/ = l[/i(i? -R)- ft'] , (2.20) 
and the off-diagonal terms vanish. The effective energy density 

Peff = P + PT ^ P + ftt, (2.21) 
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which is deduced from general relativity, p^ff 
the effective equation of state is 

rpr 



W. 



P + Ttt 

which is induced by the dynamic torsion. 



Pt = is an effective pressure, and 



(2.22) 
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Figure 1. Evolution of the Hubble parameter H, the temporal component of the 
torsion $, the affine scalar curvature R and the deceleration parameter q as functions 
of time. We have chosen the parameters ai = 0.07, b = 1.50 and the initial values 
H{1) = 1, $(1) = 1.50, R{1) = 1.20. 



In the case of Oi > 0, Nester et al showed that the scalar mode can contribute 
an oscillating aspect to the expansion rate of the universe [7]. This oscillatory nature 
can be illustrated in Fig. [T] where we have chosen oi = 0.07, b = 1.50, H{1) = 1, 
$(1) = 1.50, -R(l) = 1.20 and set the current time to = 1- According to scaling fl2.18p . 
the present value of the Hubble parameter is unity. Obviously, H is damped-oscillating 
at late time and q = —0.33 is negative today, which means the expansion of the universe 
is currently accelerating. The value of q turns from positive to negative when the time 
is around t ^ 0.51, which is the epoch the universe began to accelerate. 

However, the above result is dependent on the choice of initial data and the values 
of the parameters. Then, the scalar torsion cosmology is unsuited to solving the fine- 
tuning problem in the case of ai > 0. In the following sections, we'll investigate the 
statefinder and give the dynamics analysis for all ranges of the parameter ai. 



3. Statefinder diagnostic 

For the spatially fiat ACDM model the statefinder parameters correspond to a fixed 
point {1, 0} while {1, 1} for the standard cold dark matter model (SCDM) containing no 
radiation. Since the torsion cosmology have used the dynamic scalar torsion (a geometry 
quantity in the Riemann-Cartan spacetime), the torsion accelerating mechanism is 
bound to exhibit an essential distinction in contrast with various dark energy models. 
Therefore, its statefinder diagnostic is sure to reveal differential feature. Let us now 
study the torsion cosmological model in detail. Using Eqs. fl2.14p - fl2.17l) . we have the 
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deceleration parameter 




(3.1) 



and the statefinder parameters 

r = 1 + (6/i + bR) \bH{36H^ -3R + 4$^ - 24/7) - 12$] {lOSaibH^)' 



(3.2) 



and 



s = 



4(6/i + bR) [4bH{3H - $)2 - 3bHR - 12^$] 



(3.3) 



3bH [(6/i + bR){36H^ - 2AH^ + 4<l>2 - 3R) - 54fiR] ' 



In the following we will discuss the statefinder for four differential ranges of the 
torsion parameter ai, respectively. Firstly, we consider the time evolution of the 
statefinder pairs {r,s} and {r,q} in the case of ai > 0. In Fig. [21 we plot evolution 
trajectories in the q — r and s — r planes, where we have chosen ai = 1/2 and 6 = 4. 
We see easily that cosmic expansion alternates between deceleration and acceleration in 
the evolving trajectories of g — r plane, and the amplitude becomes larger and larger as 
increase of time. The trajectories in the s — r plane is quite complicated, so we mark its 
sequence by the ordinal number. Every odd number curve evolves from finite to infinite, 
but even number curve evolves from infinite to finite. These are quasi-periodic behaviors 
which corresponds to the numerical solution of Ref. [7]. Noticeably, the trajectories will 
never pass ACDM point {1,0}. 

Secondly, we discuss the time evolution of the statefinder pairs {r, s} and {r, q} for 
the case of — | < ai < 0. We plot evolving trajectories in Fig. [3], where we have chosen 
O'l = —jq and b = 3. We see clearly that the cosmic acceleration is guaranteed by the 
dynamic scalar torsion in the evolving trajectories of g — r plane, and the curves will 
converge into ACDM point. The evolving trajectories go through a climbing-up stage 
first, then get into a rolling-down stage in the s — r plane. Lastly, trajectories tend to 
ACDM point {1,0}. Furthermore, the only one forms a loop that starts from {1,0} 
then evolves back to {1,0}, and others show in the shape of tadpole. 
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Figure 2. Evolving trajectories of the statefinder in the q~r and s — r planes for the 
case of torsion parameter ai > 0, where we choose the parameters oi = ^ and 6 = 4. 
The arrows show the direction of the time evolution. 
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Figure 3. Evolving trajectories of the statefinder in the q — r and s — r planes for 
the case of torsion parameter — ^ < ai < 0, where we choose the parameters ai ^ —-^ 
and 6 = 3. The arrows show the direction of the time evolution. 




-0 25 -0.2 -0.1,5 -0.1 -0.05 -6 -5 -4 -3 -2 -1 

1/q 1/s 



Figure 4. Evolving trajectories of the statefinder parameters in the 1/q — 1/r and 
1/s — 1/r planes for the case of —1 < cii < — i, where we choose Oi = — ^ and 6 = 8. 
The arrows show the direction of the time evolution. 
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Figure 5. Evolving trajectories of the statefinder in the q — r and s — r planes for 
the case of torsion parameter ai < —1, where we choose the parameters ai = —2 and 
6 = 4. The arrows show the direction of the time evolution. 

Thirdly, we discuss the time evolution of the trajectories for the case of — 1 < ai < 
— |. We plot evolving trajectories in Fig. HI where we have chosen ai = —j and b = 8. 
Obviously, the cosmic acceleration can happen since deceleration parameter is negative. 
|g|, |s| and r become larger and larger first, then less and less as the cosmic time increase. 

Finally, we consider the time evolution of the statefinder pairs {r, s} and {r, q} in 
the case of ai < —1. In Fig. [5l we plot evolving trajectories in the q — r and s — r planes. 
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where we have chosen ai = —2 and h = A. We find easily that the evolving trajectories 
analogous to the case of ai > except trajectories pass the ACDM point. 

To sum up, it is very interesting to see that the scalar torsion naturally provide 
the accelerating force in the universe for any torsion parameter Oi. However, it is 
dependent on torsion parameters that there is a decelerating (g > 0) expansion before 
an accelerating (g < 0) expansion. The statefinder diagnostics show that the universe 
naturally has an accelerating expansion at low redshifts (late time) and a decelerating 
expansion at high redshifts (early time) for the cases of — | < ai and ai < —1. Obviously, 
scalar torsion cosmology can avoid some of the problems which occur in other models. 
If we refuse the possibility of non-positivity of the kinetic energy, we will employ normal 
assumption, i. e., ai > 0. In this case, the effect of torsion can make the expansion 
rate oscillate. With suitable adjustments of the torsion parameters, it is possible to 
change the quasi-period of the expansion rate as well as its amplitudes. It is worth 
noting that the true values of the statefinder parameters of the universe should be 
determined in model-independent way. In principle, vq, sq can be extracted from some 
future astronomical observations, especially the SNAP-type experiment. 

Why there are new features for the statefinder diagnostic of torsion cosmology? 
Why the torsion parameter ai is divided into differential ranges by the statefinder 
diagnostic? Answer is very simple. In fact, the statefinder diagnostic has a direct bearing 
on the attractor of cosmological dynamics. Therefore, we will discuss the dynamic 
analysis in next section. 

4. Dynamics analysis 

Eqs. f l2.14p -( !2.16p is an autonomous system, so we can use the qualitative method of 
ordinary differential equations. Critical points are always exact constant solutions in the 
context of autonomous dynamical systems. These points are often the extreme points of 
the orbits and therefore describe the asymptotic behavior. If the solutions interpolate 
between critical points they can be divided into a heteroclinic orbit and a homoclinic 
orbit (a closed loop). The heteroclinic orbit connects two different critical points and 
homoclinic orbit is an orbit connecting a critical point to itself. In the dynamical analysis 
of cosmology, the heteroclinic orbit is more interesting [19] . If the numerical calculation 
is associated with the critical points, then we will find all kinds of heteroclinic orbits. 

According to equations fl2.14p - fl2.16l) . we can obtain the critical points and study 
the stability of these points. Substituting linear perturbations R = Rc+SR, $ = $c+(^^ 
and H = Hc + 6H near the critical points into three independent equations, to the first 
orders in the perturbations, gives the evolution of the linear perturbations, from which 
we could yield three eigenvalues. Stability requires the real part of all eigenvalues to be 
negative. There are five critical points {He, $c; Rc) of the system as follows 



(i) (0,0,0) 




3(1 + ai) 
8 




6(1 + ai) 
b 



) 
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where A = V'ii^+'^W*"»\ g _ _ (i+a.)(5+»a.) C = _ V°Ki+°.)'(i+a-) . 

8ai(l+ai) ' aio a^o 

corresponding eigenvalues of the critical points (i)-(v) are 



(4.1) 
The 



(i)(0, 



:ii) -v5 + c,-3 



2(1 + ai) / 2(1 + ai) . 
aih ' V o-ih 
■3(1 + 01) 



111 



IV 



'v/5TC, 3 (^(1±^ - a) V^TC, - 3a) v/bTc) 



'3(1 + ai) \ r- /l + 9ai 



-3 ^ ' ' + A VW^, - ' + 3A 



(V) (v/B^, 3 (^(1±^ + a) v/B^c, + 3a) v/B^c) (4.2) 

Using Eq. (14.11) . we find that there is only a critical point (0,0,0) in the case of 
ai > 0. From Eq. (14.21) . the corresponding eigenvalue is (0, —\f^i-, \f^i)i so (0, 0, 0) 
is an asymptotically stable focus. If we consider the linearized equations, then Eqs. 
are reduced to 

H = ^R. ^ = i^R. R = -^^- (4-3) 

DOi 2ai 

The linearized system (14.31) has an exact periodic solution 

H = -aRo sin cut + -$o cos oot + Hq - -^0, 
3 3 

$ = -(3~^Ro sin iut + $0 cos ut 

R = Ro cos ut + /3$o sin ut (4.4) 

where u = -y^, a = /5 = and ifo = H{0), % = m and Ro = R{0) 

are initial values. Obviously, [H, 0,0) is a critical line of center for the linearized Eqs. 
(14. 4p . In other words, there are only exact periodic solutions for the linearized system, 
but there are quasi-periodic solutions near the focus for the coupled nonlinear equations. 
This property of quasi-periodic also appears in the statefinder diagnostic with the case 
of ai > 0. 

In the case of —1/9 < ai < 0, the critical point (ii) is a late time de Sitter 
attractor and (iii) is a negative attractor. The properties of the critical points are 
shown in table [H The de Sitter attractor indicates that torsion cosmology is an elegant 
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Table 1. The physical properties of critical points for — ^ < < 



critical points property ^e// stability 



(i) saddle — oo unstable 

(ii) positive attractor -1 stable 

(iii) negative attractor -1 unstable 

(iv) saddle -1 unstable 

(v) saddle -1 unstable 



Table 2. The physical properties of critical points for ai < — 1 



critical points property "^e// stability 



(i) focus ±oo stable 

(ii) saddle -1 unstable 

(iii) saddle -1 unstable 

(iv) saddle -1 unstable 

(v) saddle -1 unstable 



scheme and the scalar torsion mode is an interesting geometric quantity for physics. In 
the dynamical analysis of cosmology, the heteroclinic orbit is more interesting. Using 
numerical calculation, we plot the heteroclinic orbit connects the critical point case (iii) 
to case (ii) in Fig. El This heteroclinic orbit is just corresponding to the loop in Fig. 
[3l which is from ACDM point to ACDM point. Furthermore, the trajectories with the 
shape of tadpole correspond to saddles. 




Figure 6. The phase diagrams of (iJ, $,i?) with —1/9 < ai < 0. The heteroclinic 
orbit connects the critical points case (iii) to case (ii). We take ai = —1/10, 5 = 4. 
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Figure 7. The phase diagrams of {H,^,R) with ai < —1. We take ai = —2,6 = 4 
and tlic initial value (5, 1, ||). (0, 0, 0) is an asymptotically stable focus point. 

In the case of — 1 < ai < — |, there is only an unstable saddle (0,0,0) where the 
effective equation of state Wg// tends to — oo. Therefore, the trajectories in Fig. Hlshow 
that |g|,|s| and r become larger and larger, then less and less as time increases. 

In the case of ai < —1, the properties of the critical points are shown in table [2j 
The trajectories correspond to the stable focus (see Fig. 6) and unstable saddles with 
Weff = —1. Therefore, the trajectories pass through the ACDM point. 

5. Fit the torsion parameters to supernovae date 

In Ref. [7] , the authors have compared the numerical values of the torsion model with 
the observational data, in which they fixed the initial values Hq, $o and Rq, and torsion 
parameters ai and b. In this section, we fixed the initial value, then fit the torsion 
parameters to current type la supernovae data. 

The scalar torsion cosmology predict a specific form of the Hubble parameter H{z) 
as a function of redshifts z in terms of two parameters ai and b when we chose initial 
values. Using the relation between diiz) and the comoving distance r{z) (where z is 
the redshift of light emission) 

dL{z) = r{z){l + z) (5.1) 

and the light ray geodesic equation in a fiat universe cdt = a{z)dr{z) where a{z) is the 
scale factor. 

In general, the approach towards determining the expansion history H{z) is to 
assume an arbitrary ansatz for H{z) which is not necessarily physically motivated but is 
specially designed to give a good fit to the data for di^z). Given a particular cosmo logical 
model for Hi Z] Ci\ , . . . , dn ) where ai, a„ are model parameters, the maximum likelihood 
technique can be used to determine the best fit values of parameters as well as the 
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goodness of the fit of the model to the data. The technique can be summarized as follows: 
The observational data consist of apparent magnitudes mi{zi) and redshifts Zi with 
their corresponding errors Um. and cr^.. These errors are assumed to be gaussian and 
uncorrelated. Each apparent magnitude mj is related to the corresponding luminosity 
distance cii by 

dL{z) 



m{z) = M + 5 log 



10 



Mpc 



+ 25, 



(5.2) 



where M is the absolute magnitude. For the distant SNela, one can directly observe 
their apparent magnitude m and redshift z, because the absolute magnitude M of 
them is assumed to be constant, i.e., the supernovae are standard candles. Obviously, 
the luminosity distance d^lz) is the 'meeting point' between the observed apparent 
magnitude m{z) and the theoretical prediction H{z). Usually, one define distance 
modulus z/(z) = m{z) — M and express it in terms of the dimensionless 'Hubble-constant 
free' luminosity distance defined by -0^(2) = HQdL{z)/c as 

iy{z) = 5 \og,,{DL{z)) + z/o, (5.3) 

where the zero offset uq depends on Hq (or h) as 



1/0 = 5 log 



10 



Mpc 



25 



51ogio/i + 42.38. 



The theoretically predicted value D*^{z) in the context of a given model H{z] ai, 
can be described by [20] 

Ho 



Dfiz) 



l + z) 



dz' 



(5.4) 
..,an) 

(5.5) 



H{z'; ai, ...an) 

Therefore, the best fit values for the torsion parameters (ai,6) of the model are found 
by minimizing the quantity 

2 



N 



X^(ai,6) = J2 



[u'^'iz,] - 5\og,oDf{zf,a^,b) - 1^0 



i=l 



(5.6) 



Since the nuisance parameter z/q is mo del- independent, its value from a specific good fit 
can be used as consistency test of the data [21] and one can choose a priori value of 
it (equivalently, the value of dimensionless Hubble parameter h) or marginalize over it 
thus obtaining 



X^(ai,6) = A{ai,h) 



Bjai.bf 
C 



In 



2^ 



where 



A{aub) = 



N ^jy 
i=l 



[u'''^{z,)-5\og,,Df{zf,a^,b) 



N 



B{aub)=Y: 



[u"'>'iz,)-5\og,,Df{z,;ai,b) 



(5.7) 

(5.8) 
(5.9) 



1=1 
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and 



N 1 



c = E-2- (5-10) 



i=l 



In the latter approach, instead of minimizing b), one can minimize x^{ai, b) which 

is independent of uq. 

The Eqs. f l2.8ti2.10]) can be solved explicitly by a series in the form 

oo 

H{t,to) = Ho[l + Y.hnu''], (5.11) 

n=l 

where u{t) = e'^^^~^°'^ — 1 and 

Hq Rq i?o bHoRo bRl 



hi 



ijj 2ujai GujHq GcuaiHo 12uai lAAuaiHi 







^ $0 ^ ^ ^ bRo% _ $g _ _ bRo<l>l 







h ^ ^0 I Hp 2hiHo hiHp bhiRpHp ^ bRl ^ hi% 
UJ Aujoi UJ 2ujai 12ujai 288ujaiHQ 6uj 
hi<i>o , bRo% , bhiRo% '^l 6i?o<fo 



+ 



Gojai 36u;ai 36u;ai 36ujHq 36ujaiHo lOSujaiHo ' 



Using the general relation between Hubble parameter H{t) and the redshift z, z can be 
written function of t 

z{t) = exp ij^^' H{t)dt\ - 1, (5.13) 

However, the convergence radius of the series (15. lip is | ^ | , so we can use the expansion 
directly in the case of the redshift being z < z{tQ + '^). By the numerical calculation, 
we find that t* ~ to + ^ corresponds to ~ 0.45 for the valuses of parameters Rq and 
$0 in the Fig. [9l For the case of 2; > 2*, we should use a direct analytic continuation. 
Weierstrass [22] had built the whole theory of analytic functions from the concept of 
power series. Given a point ti = + at* (0 < a < 1), the function H{t) has a Taylor 
expansion 

n 

H{t,ti) = Hiti,to)[l + E /i„(e-(*-*^) - 1)"]. (5.14) 

n=l 

where the coefficients hn is still expressed as Eq. f l5.12p and Hiti.to) can be decided 
by Eq. (15. lip . The new series defines an analytic function H(t,ti) which is said to be 
obtained from H{t, to) by direct analytic continuation. This process can be repeated 
any number of times. In the general case we have to consider a succession of power 
series if(t,to), H{t,ti),...,H{t,tjn), each of which is a direct analytic continuation of the 
preceding one. By using this method we have the evolution of Hubble parameter H{t). 
Furthermore, we have the function H{z) from Eq. (I5.13p . In fact, we need only to 
consider the case of 2; < 2 for the ESSENCE supernovae data. 
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We now apply the above described maximum likelihood method using the 
ESSENCE supernovae data which is one of the reliable published data set consisting of 
192 SNela (A^ = 192). Beside the 162 data points given in table 9 of Ref. [23], which 
contains 60 ESSENCE SNela, 57 SNLS SNela and 45 nearby SNela, we add 30 SNela 
detected at 0.216 < z < 1.755 by the Hubble Space Telescope [21] as in Ref. [25] . 

In table [3|, we show the best fit of torsion parameters at different initial values of Rq 
and $0- In Fig- El contours with 68.3%, 95.4% and 99.7% confidence level are plotted, 



Table 3. The best fit of torsion parameters for different initial values of Rq and $o 


Rq 




ai 


b 


0.25 


0.35 


-0.10 


1.44 


0.20 


0.34 


-0.08 


1.80 


0.15 


0.34 


-0.06 


2.40 


0.10 


0.33 


-0.04 


3.60 



in which we take a marginalization over the model-independent parameter z/q. The best 
fit as showed in the figure corresponds to Oi = —0.06 and b = 2.40, and the minimum 
value of = 355.68. For ACDM, one can get x'^i^CDM) = 355.74 and the best 
fit Qmo = 0.26. Therefore, it's easy to know that ACDM is consistent at the la level 
with the best fits of scalar torsion cosmology. In Fig. [HI we show a comparison of the 
ESSENCE supernovae data along with the theoretically predicted curves in the context 
of scalar torsion and ACDM . We can see that the scalar torsion model(-Ro = 0.15, 
$0 = 0.34, ai = —0.06, b = 2.4) gives a close curve behavior to the one from ACDM 
{Qmo = 0.26, f^A = 0.74). Clearly, the allowed ranges of the parameters ai and b favor 
the case of — ^ < ai < if we chose Rq = 0.15 and $o = 0.34. 




0.5 1 1.5 2 1.2 1.4 1.6 1.8 2 2.2 



Figure 8. Comparison of scalar torsion model and ACDM with the ESSENCE 
supernovae data via the relation between the redshift z and the distance modulus 
/i. The scalar torsion model(i?o — 0.15, $o = 0.34, ai ~ —0.06, b — 2.4) gives a close 
curve behavior to the one from ACDM {^Imo — 0.26, Ha = 0.74). 
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Figure 9. The 68.3%, 95.4% and 99.7% confidence contours of torsion parameters ai 
and b using the ESSENCE SNela dataset. Here we have assumed Hq = I, Rq = 0.15 
and $0 = 0.34. 

6. Conclusion and discussion 

We have studied the statefinder diagnostic to the torsion cosmology, in which an 
accounting for the accelerated universe is considered in term of a Riemann-Cartan 
geometry: dynamic scalar torsion. We have shown that statefinder diagnostic has a 
direct bearing on the critical points. The statefinder diagnostic divides the torsion 
parameter ai into four ranges, which is in keeping with the requirement of dynamical 
analysis. Therefore, the statefinder diagnostic can be used to an exceedingly general 
category of models including several for which the notion of equation of state is not 
directly applicable. The statefinder diagnostic has the advantage over the dynamical 
analysis at the simplicity, but the latter can provide more information. 

The most interesting characteristic of the trajectories is that there is a loop in the 
case of — I < ai < 0. This behavior corresponds to the heteroclinic orbit connecting the 
negative attractor and de Sitter attractor. The trajectories with the shape of tadpole 
show that they pass through the ACDM fixed point along the time evolution, then the 
statefinder pairs are going along with a loop and they will pass through the ACDM fixed 
point again in the future. It is worth noting that there exists closed loop in the Ref. 
[13], but there is no closed loop which contains the ACDM fixed point. These behaviors 
indicate that torsion cosmology is an elegant scheme and the scalar torsion mode is an 
interesting geometric quantity for physics. Furthermore, the quasi-periodic feature of 
trajectories in the cases of ai > or ai < — 1 shows that the numerical solutions in 
Ref. [7] are not periodic, but are quasi-periodic near the focus for the coupled nonlinear 
equations. 

We fixed only the initial values, then fitted the torsion parameters to current SNela 
dataset. We find that the scalar torsion naturally explain the accelerating expansion 
of the universe for any torsion parameter ai. However, it is dependent ai and b that 
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there is a decelerating expansion before an accelerating expansion. The statefinder 
diagnostics show that the universe naturally have an accelerating expansion at late time 
and a decelerating expansion at early time for the case of ai > — | and ai < — 1. 
If we refuse the possibility of non-positivity of the kinetic energy, we have to employ 
normal assumption (ai > 0). Under this assumption, the effect of torsion can make 
the expansion rate oscillate. Furthermore, with suitable adjustments of the torsion 
parameters and initial value, it is possible to change the quasi-period of expansion rate 
as well as its amplitudes. In order to have a quantitative understanding of the scalar 
torsion cosmology, the matter density p, the effective mass density pe// = P + Pt, and 
the quantity 3pT + Peff are important. This scenario bears a strong resemblance to 
the braneworld cosmology in a very different context by Sahni, Shtanov and Viznyuk 
[26]. The Qm parameters in the torsion cosmology and in the ACDM cosmology can 
nevertheless be quite different. Therefore, at high redshift, the torsion cosmology 
asymptotically expands like a matter-dominated universe with the value of inferred 
from the observations of the local matter density. At low redshift, the torsion model 
behaves like ACDM but with a renormalized value of Q^^'^'^ . The difference between 
flm and Q^^^^ is dependent on the present value of statefinder parameters ro, sq. A 
more detailed estimate, however, lies beyond the scope of the present paper, and we 
will study it in a future work. Finally, tq and Sq should be extracted from some future 
astronomical observations in principle, especially the SNAP-type experiments. 
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